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Abstract 
We show that self-annihilating neutralino WIMP dark matter accreted onto neutron stars may provide a mechanism to seed 
compact objects with long-lived lumps of strange quark matter, or strangelets, for WIMP masses above a few GeV. This effect 
may trigger a conversion of most of the star into a strange star. We use an energy estimate for the long-lived strangelet based 
on the Fermi gas model combined with the MIT bag model to set a new limit on the possible values of the WIMP mass that 
can be especially relevant for subdominant species of massive neutralinos. 


There is compelling evidence [1] that most of the mat- 
ter in the Universe is dark matter (DM). Amongst the 
beyond-the-Standard Model particles that may consti- 
tute DM, WIMPs (weakly interacting massive particles), 
typically the LSP neutralinos of SUSY models, are the 
most frequently considered. WIMPs from the galactic 
halo can be accreted onto compact massive objects |2], 
such as neutron stars (NS) ia, and white dwarfs ap 
Our focus here will be on the consequences of neutralino 
WIMP capture by NS. Accretion of DM has been shown 

,\q] to play a role in late cooling for NS over time-scales 
longer than Tys > 107 yr when the WIMP annihilation 
rate has already reached equilibrium with the accretion 
rate. Then the rate of released energy becomes time and 
temperature-independent, and will dominate over photon 
cooling processes. However this effect is hard to detect. 

The signature of WIMP self-annihilations (assuming 
WIMPs are Majorana particles) is energy release in the 
range equivalent to twice the WIMP mass, leading to sub- 
sequent particle production and/or heating of the star. 
Current predictions of the WIMP mass span the range 
from 1 GeV/c? up to 10 TeV/c? [S]. Other, more exotic, 
DM candidates include super-heavy WIMPzillas |g, 
with masses in the range 101° — 10!° GeV/c? and SIM- 
Pzillas - heavy and strongly interacting 12) (see re- 
cent experimental constraints from ICECUBE |13}). 

Compact objects such as black holes or NS are consid- 
ered to be gravitational accretors of DM. Non-rotational 
NS models, yielding gravitational mass, M, and radius, 
R, are based on the Tolman-Oppenheimer-Volkoff (TOV) 
equations, given an equation of state (EoS) that relates 
pressure to energy density and temperature, p = p(e,T). 
A typical NS has a mass M% 1.4 Mo, radius R ~ 10 km, 
surface temperature Tsuf œ% 100 keV and central mass 
density p © 1.5x 10!° g/cm. There are a variety of 
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model predictions for NS composition [14] but there is 
as yet no consensus about the fundamental nature of the 
matter in NS interiors (15). 

Witten proposed in 1984 that strange quark mat- 
ter (SQM) made of u, d and s quarks is absolutely sta- 
ble and forms the true ground state of hadronic matter, 
and that strange stars, made of SQM and bound by the 
strong interaction, should exist with very different prop- 
erties from those predicted for hadronic NS aj. The 
mechanism of conversion of metastable, charged neutral, 
neutron star matter in beta equilibrium, formed due to 
the leptonic weak interaction p + e © n +v, to a two- 
flavor quark matter (ud matter) and subsequently to uds 
matter is not completely understood, and a variety of 
models have been suggested (see e.g. H23). Here we 
show that WIMPs may trigger a conversion of hadronic 
NS matter to SQM through an external mechanism of 
seeding the NS with long-lived lumps of SQM as the re- 
sult of energy release from DM self-annihilations in a dis- 
tribution of gravitationally-captured WIMPs in the NS 
core. 

The accretion rate of WIMPs captured by a typical NS 
is given by (2), a 


3.042 x 10% PDM ; 1 


F ae EEn Á 
my(GeV) PDM,0 


(1) 
assuming a WIMP-nucleon interaction cross section 
oyn > 10745 cm?. We assume that the local DM mass 
density in the vicinity of the NS in an average galaxy 
is similar to that in the solar neighbourhood, about 
ppM,o = 0.3 GeV/cm?, and that all incoming WIMPs un- 
dergo one or more scatterings while inside the star. Here 
my is the mass of the WIMP. Due to competing effects 
of annihilation and evaporation, the number of accreted 
WIMPS at time t is obtained by solving the differential 
equation 


N =F- Tannih — Pevap; (2) 


where Tannin is the self-annihilation rate calculated as 
Tannih = (JannihV) J naV = Ca N(t)?, (Jannih) is the 


product of thermally-averaged WIMP self-annihilation 
cross section and velocity, and ny = at the number 
density of WIMPs inside the NS, is assumed to be con- 
stant. The evaporation rate, evap, decays exponentially 
with temperature ~ e~GM™x/RT | and is negligible with 
respect to the annihilation rate for NS with internal T ~ 
0.1 MeV [25]. With this simplification, the population of 
WIMPs at time t is given by 


N(t) = (Fr)tanh(t/r), (3) 


where the time-scale is r = 1/,/FCa. For t >> 7, when 
the equilibrium between accretion and annihilation has 
been reached, the number of particles accreted is time- 
independent, N = Fr. 

Assuming a regime when velocities and positions of the 
WIMPs follow a Maxwell-Boltzmann distribution with 
respect to the centre of the NS, the thermalisation volume 


in the compact star has a radius |5] rin = ( To \i/2 = 
14 3 2 
64 (api )1/2( BL 1/2 Oe 1/2 (cm). Taking 


typical NS conditions as central internal temperature 
Te = 10° K, p-=10!° g/cm? and my = 1GeV/c? — 
10 TeV /c?, tin © 2 x 10° — 2 x 10?cm. Then the ther- 
malization volume is Vin = anren x (107? — 10-8) Vys. 

The exact energy released due to WIMP. self- 
annihilation is dependent on the nature of WIMPs and 
on the output product channel and has to be calculated 
separately for each case. As we are interested in order 
of magnitude estimates, we assume that after an anni- 
hilation occurs, a fraction of this energy, given by the 
efficiency factor le, f, is deposited locally. The rate of 
energy released in annihilation processes in the star in 
the thermal equilibrium volume can be expressed as 


Eannih = CaN? m,c? = JFmg®, (4) 


where f ~ 0.01 — 1. If most of the total energy is not 
converted into neutrinos, then f is close to unity. 

Eq. (4) includes the option that ppm, in the vicinity of 
the NS, may differ from the local standard DM density, 
ppm,o. Taking ppw=ppmM,o and an efficiency rate f = 
0.9, we have that in the range my ~ 1 — 104 GeV/c? 


Eannih = 2.74 x 107° — 107? GeV/s. (5) 


The energy released will convert partly into heat, which 
can stimulate u,d quark bubble formation by thermal 
fluctuations in metastable hadronic matter via strong in- 
teractions (20). The timing and conditions of these tran- 
sitions depend strongly on the pressure in the central 
regions of the star given by the EoS of hadronic mat- 
ter. Then ud matter undergoes non-leptonic weak reac- 
tions, such as u + d + u + s, to form drops of strange 
uds matter, which has lower energy as a result of the 
reduction in Fermi energies through the introduction of 
a new flavor. The stability, among other properties of 


these drops, commonly called strangelets, has been ex- 
tensively studied (see e.g. l) and depends on the elec- 
trical charge, strangeness fraction and size. The energy 
needed to form a stable long-lived strangelet of baryonic 
number A can be calculated from its quark constituents, 
in a first approximation, using either the MIT bag model 
with shell mode filling or the liquid drop model (for de- 
tails see e.g. Bd). The former approach is important 
particularly at low A but becomes impractical for large 
A. Minimum values of A and lifetimes of strangelets 
have been calculated in the MIT bag model approach 
at zero temperature for different values of the bag 
constant, B, charge and strangeness. Finite tempera- 
ture effects on strangelets have also been examined [30}. 
Although they do not appreciably modify the minimum 
baryon number of a strangelet, they shift stability condi- 
tions towards the high negative electrical charge fraction 
and high strangeness. 

In the ideal Fermi-gas approximation, the binding en- 
ergy of a strangelet with baryonic number A composed 
of massive quarks of flavor i (i = u, d, s) is 


EA (mi, mi, B) = XO (0 + Napi) + BV (6) 


i 


where N; is the quark number in the strangelet of baryon 
number A = 3(Nu + Na+ Ns), pilna) is the quark chem- 
ical potential at baryonic number density na = DD Nij. 


ni are the ith-quark number densities in the strangelet. 
The thermodynamical potential for the ith-quark type 
with mass m; is given by 0; (4i,m:) = OQ VV + 24,95 + 
QicC being V = A/na = sn R3 the volume, S = 47R? 
the surface and C = 87R the curvature of the spheri- 
cal strangelet. Expressions for these potentials can be 
obtained from . The masses of the quarks are taken 
as My, = 2.55 MeV, ma = 5.04 MeV, m, = 104 MeV 
respectively igh. 

In order to consider charged strangelets a Coulomb 


4 Z2 2 
correction term can be added as Ecoul = i(i + sZ) 


with Zy = 5 qin; V and Z = 5 qi given the ith-quark 


charge qi. Then, the energy of a strangelet with baryon 
number A is 

Eiet = EA (mi, mi, B) + Ecoul- (7) 
Note that we have neglected a zero-point motion correc- 
tion as done in |28]. Estimation of the minimum value of 
A in a long-lived strangelet is model-dependent but typ- 
ical minimum A values are in the range Amin œ 10 — 600 
BJ. Schaffner et al. find strangelets with ‘magic’ 
numbers of quarks for bag constant values 145 < B 1/4 < 
170 MeV. If smaller strangelets than the minimum A are 
created, they will decay rapidly. However, it remains for 
further investigation to see whether quark drop percola- 
tion plays a significant role here. This may happen as a 
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FIG. 1: Mass of the DM WIMP as a function of B'/* for dif- 
ferent values of the baryonic density, na, taking the efficiency 
factor f = 0.9. 


result of multiple WIMP annihilations over scales of sev- 
eral fm (typical size of a strangelet). A long-lived cluster 
may be formed from different small clusters with interdis- 
tance l on the diffusion time-scale, tai ~ l/vaig¢, before 
they decay over a mean lifetime, r ~ 10719 — 10~° s. 
This effect would enhance the seeding scenario presented 
here by increasing the energy available for strangelet for- 
mation. The more elaborate calculations show that 
among all possible metastable strangelets, the long-lived 
ones have lifetimes of days. Since this time-scale is larger 
than the conversion time-scale teony % 100s [21], it is in 
principle possible that strange stars could be formed if 
this conversion is triggered. As an estimate, the rate 
of formation of long-lived strangelets with, for example, 
A = 10 would be 


Niiet = Pah Efa- (8) 


For f = 0.9, na = nao, My = 1 GeV/c?, a number 
Netet © 1028 s7! are created. Assuming that in the cen- 
ter of the star, self-annihilations are dominant, the en- 
ergy deposited in the medium, per annihilation, will be 
2fm,c?. This allows us to set a limit on the energy scale, 
at a given central baryonic density na. 


2fmyc? 2 Exe (ui(na), Mi, B). (9) 


Note that the mass of the WIMP can be related to the 
strong interaction through the MIT bag constant. Eq. (0) 
provides a strong limit to complement current estimates 
from direct as well as indirect searches. In Fig[i] we 
show the WIMP mass as function of B'/* for a long-lived 
strangelet baryon number A = 10 and an efficiency rate 
f = 0.9 for na = nao (solid line), na = 2n4,o (long 
dashed line) and na = 5n4,o (short dashed line), with 
nA, o = 0.17fm7?. It can be seen that there is a weak 
dependence on nuclear physics input, parameterized by 
B. As the central baryon density increases, the lower 
limit on the WIMP mass decreases below ~ 4GeV/c?. 
We adopt na = 5n 4,0 as the highest value of the central 
baryon density. 
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FIG. 2: WIMP mass as a function of long-lived strangelet 
baryon number A for (a) limiting values of B'/* = 140,170 
MeV, and na/na,o = 1,3,5 (b) BY/4 = 170 MeV, 3na,o and 
f = 0.1, 0.6, 0.9. 


In Fig. [2] we show (a) the WIMP mass as a func- 
tion of A in the long-lived strangelet for an efficiency 
rate f = 0.9 and different values of the central den- 
sity na = nao (dash-dotted line), 3n4,9 (thick solid 
line), 5n4,0 (thin solid line) for limiting values of the 
bag constant. For each pair of curves, the lower line is 
for B!/4 = 145 MeV and the upper line is for B!/4 = 170 
MeV. For a typical baryonic density of ~ 3n4,0, a lower 
linit for the WIMP mass my Z 4 GeV/c? is predicted. 
In Fig[2] we show in (b) the WIMP mass as a function of 
the long-lived strangelet baryon number A in the range 
10— 16 according to ref. [29] for different values of the effi- 
ciency rate f and B*/4 = 145 MeV and na = 3nao. We 
consider f = 0.1 (thin solid line), f = 0.6 (dashed line) 
and f = 0.9 (thick solid line). As more massive WIMPs 
annihilate, the smaller the efficiency rate becomes that is 
required to trigger a conversion from nucleon to strange 
quark matter. 


In summary, through the mechanism of self- 
annihilation of DM candidates in the central regions of 
a typical NS proposed in this work, we have derived a 
limit on the mass of DM particle candidate using the 
fact that there is a minimum long-lived strangelet mass 
needed to trigger a conversion from nuclear matter into 
strange quark matter. In the range my Z 4 GeV/c? the 
energy released in WIMP self-annihilation is sufficient to 
burn nucleon matter into a long-lived strangelet that will 
trigger full conversion to a strange star. It is important 


to emphasize that so far the lack of knowledge of DM 
properties prevents us from making any definite state- 
ment on the WIMP mass other than excluded regions, 
but this work presents a scenario compatible with current 
experimental direct and indirect searches. Two impor- 
tant pieces of observational information would provide 
supporting evidence for the mechanism of NS conversion 
proposed here and potentially improve the derived limit 
on the WIMP candidate: (i) observation of a quark star 
and measurement of its mass and radius, which could 
yield constraints on its EoS; (ii) identification of forma- 
tion of strangelets and its properties at LHC or RHIC. 

In the astrophysical scenario, one consequence of the 
NS conversion process could be related to the emission 
of y-ray bursts (GRB). According to Ma et al. 
the energy release in a Super-Giant-Glitch due to the 
change in radius of an initial to final star configuration 
is AE ~ 10°342 erg. The calculated rate of GRBs in 
the galaxy is Rapp ~ 107% yr~! in good agreement with 
observations. Such a huge burst would happen only once 
in the lifetime of the NS. Bombaci et al. , pointed out 
that the energy liberated during the NS conversion pro- 
cess will be carried out mainly by neutrinos that will be 
visible near the surface as GRBs. In addition, the strong 
magnetic field in the NS can cause collimated GRBs and 
the rotation of the NS could be a source of anisotropic 
emission. They have also investigated the effect of un- 
certainties introduced by different EoS. Recent measure- 
ments of the source of short GRB 090510 by Fermi LAT 
report an energy FE = (1.08 + 0.06) x 105° erg in the (10 
keV-30 GeV) range and a time duration of most of 
the signal, T9) = 0.6 — 9 s, compatible with the expected 
energy release and duration of the conversion of a 
hadronic to quark star (21h. 

An interesting open question is whether conventional 
pulsar glitches observed so far only in hadronic stars 
would be possible in quark stars which are predicted to 
have no crust or very thin one. Some models sug- 
gest that it is unlikely that quark stars would give rise to 
glitches as these require a crust with charged strangelets 
inmersed both in electron gas and a superfluid and there 
is no superfluid in the crust of quark stars. However if 
the quark matter in the core of the star were to exist in 
crystalline color-superconducting phase, glitches might 
origite there. 

The rate of conversions of hadronic to strange NS via 
the mechanism presented in this work is directly related 
to the DM density, since it depends on the galactic DM 
profile. Thus this rate could be potentially observable as 
a function of position of the NS in the galaxy. If WIMP 
annihilation is considered as an effective way of seeding 
NS with strangelets, multicomponent DM may have ob- 
servable consequences in the conversion of NS to strange 
stars even if the annihilating component (e.g. SIMP2zil- 
las) is subdominant, as shown by the ICECUBE con- 
straint. It would be also interesting to explore whether 


the annihilation of DM could cause changes in superfluid 
phases of matter inside NS and its effect on cooling and 
rotation patterns. 
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